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1. 


INTRODUCTION 


Above  100  km  the  infrared  (IR)  radiation  from  2  -  50p  emitted  from 
rocket  and  missile  plumes  may  be  produced  largely  from  high  velocity  molecular 
collisions  between  IR-active  exhaust-exhaust  and  exhaust-atmospheric  species. 

At  higher  altitudes  these  TR-signatures  are  dominated  by  the  Latter  process, 
the  predominant  atmospheric  species  being  atomic  oxygen.  Since  these  inter¬ 
actions  occur  with  high  relative  energies,  collisions  between  atomic  oxygen 
and  exhaust  species  such  as  HF,  HO,  and  C09  may  result  in  excitation,  reaction, 
or  dissociation  processes.  In  the  present  work  wo  will  he  concerned  only  with 
nonreactive  excitation  processes. 

At  present.  Laboratory  data  exists  only  for  vibrational  excitation  of 
f  1  9 1  (  1 1 

C09  by  N0,  O2  and  Ar  and  H^O  by  N0 .  There  has  been  one  field  measure¬ 

ment  of  a  high  altitude  infrared  exhaust  plume  signature  from  which  an  excita¬ 
tion  cross  section  for  H.,0  and  CO.,  (by  0-atoms)  has  been  extracted^ Due  to 
the  complexity,  the  only  theoretical  work^'*  on  the  vibrational  excitation  o! 

CO^  by  0  has  considered  CO,,  as  a  diatomic  molecule. 

The  present  work  represents  an  initial  study  of  the  theoretical  characteri¬ 
zation  of  excitation  cross  sections  of  H^O  and  C0o  in  collisions  with  a 
structureless  particle,  i.e., 

^^Subbarao,  R.,  Fenn,  J.B.,  and  Kolb,  0.F,.,  .1 .  Ghent.  Phvs.  to  be  submitted. 

(2) 

Rahbee,  A.,  Gibson,  J.,  and  Dolan,  C.,  private  communication  (1980). 

^^Dunn,  M.G.,  Skinner,  G.T.,  and  Treamor,  G.F.,  AIM  J.  13,  801  (1979). 

(A) 

McIntyre,  A.,  Gersh,  M.K.,  Wheeler,  N.B.,  Frankel,  D.S.,  and  Klgin,  .1.13., 
private  communication  (1979). 

(5)Bass,  J.N.,  J.  Ghem.  Phvs.  60,  2913  (19/4). 


M  +  H20 


- ►  M  +  H  0+ 

and 

M  +  C02  - ►  M  +  C09* 

Although  a  complete  study  of  the  dynamics  of  these  processes  is  needed, 
it  has  been  shown  that  the  qualitative  concepts  of  the  information  theoretic 

/  rl  \ 

approach  are  quite  useful.  We  show  the  results  of  such  a  study.  We  find 
that  the  energy  dependence  of  these  cross  sections  appears  to  agree  quite 
well  with  that  from  the  available  experimental  results. 

The  organization  of  this  work  is  as  follows:  Section  2  presents  a  general 
outline  of  the  theoretical  methods  used,  and  Section  3  discusses  the  predicted 
excitation  cross  sections  and  compares  them  with  existing  data. 


Kor  reviews,  see  (a)  Levine,  R.D.  and  Bernstein,  R.B.,  Arc.  Chem.  Res.  7, 
193  (1974):  (b)  Bernstein,  R.B.  and  Levine,  R.D.,  Adv.  At.  Mol.  Plivs.  11, 
215  (1975);  (c)  Levine,  R.D.  and  Bernstein,  R.B.,  in  Modern  Theoretical 
Chemistry,  Vol.  II:  Dynamics  of  Molecular  Collisions,  ed .  by  W.H.  Miller 
(Plenum,  New  York,  1976),  p.  323. 


2.  THEORETICAL  METHOD 


The  information  theory  approach  to  dynamics  is  well  documented  in  the 


(  7 

literature'  ’  and  will  only  be  outlined  here. 


The  state-to-state  cross  section  for  reactants  in  quantum  state  n  and 

(5,6) 


products  in  the  final  state  n'  is  given  by 

.4 


o(n  -*■  n' ;  E)  =  — ( \  pt(E  -  En>  Pt(E  -  )  | Tnn ,  |2 


(1) 


gn  k 
n  n 


where  E  is  the  total  energy,  g.  and  E.  are  the  degeneracy  and  energy  of  the 
th  ^  t 

i  internal  quantum  level,  T  ,  is  the  transition  matrix  element  for  the 
‘  nn 


2  /  2 

n  ->  n'  transition,  k  =  2y(E  -  E  )/ti  with  y  as  the  reduced  mass,  and 

n  n  / 

Pt<E  -  E^)  is  the  translj 
ich  quantum  level,  i.e.. 


Pt(E  -  E^)  is  the  translational  density  of  states  for  a  species  in  the 


Pt(E  -  E.)  =  g.  At(E  -  E.) 


1/2 


(2a) 


with 


3/2.-1/2  2+3 
y  /2  n  h 


(2b) 


The  concept  of  a  "Prior"  cross  section,  labeled  o  ,  is  that  which  obtains  when 
all  quantum  transitions  are  taken  to  be  equally  probable,  i.e., 


nn 


constant  . 


(3) 


(7) 


Procaccia,  I.  and  Levine,  R.D.,  J.  Chem.  Phys.  64^,  808  (1976), 


3 


Inserting  Eqs.  (2)  and  (3)  into  (1)  yields 


o° (n  ->  n' ;  E) 


(4a) 


B  V°  -  Qnn’/ET> 


1/2 


(4b) 


Here,  B  is  a  collection  of  constant  factors,  the  initial  translational  energy 
is  E^,  =  E  -  E^,  and  the  energy  defect  of  the  transition,  Q  , ,  is  given  by 


Q 


nn' 


(E  ,  +  E  ,)  -  (E  +  E  )  . 
no  no 


(5) 


Note  that  E^  is  the  internal  energy  relative  to  the  zero  point  level  (for 
electronically  adiabatic  nonreactive  processes,  Eq  =  E  , ) .  It  is  interesting 
that  on  purely  kinematic  grounds  (i.e.,  no  dynamical  considerations),  the  prior 
cross  section  decreases  for  exoergic  (Qnnt  <  0)  processes  and  increases  for 
endoergic  (Q  ,  >  0)  processes. 

Now,  the  true  transition  cross  section  will  differ  from  the  prior  cross 
section  due  to  various  dynamical  effects  or  constraints.  It  is  the  goal  of 
information  theory  to  predict  the  "best"  form  of  the  cross  section  consistent 
with  the  observables  of  the  process.  This  is  accomplished  by  use  of  the  Max¬ 
imum  Entropy  Principle!^  which  simply  stated  means  that  a  system  proceeds  such 
that  the  entropy  achieves  its  maximal  allowed  value,  constrained  of  course,  to 
reproduce  the  independent  observables. 

Now,  the  state-to-state  cross  section  may  be  written  in  the  form 


o(n  v  n')  =  o(n)  P(n'|n) 


(6) 


where  o(n)  is  the  total  cross  section  out  of  state  n  and  P(n'|n)  is  the  condi¬ 
tional  probability  of  scattering  into  state  n'  given  a  reactant  in  state  n. 
namely 


4 


A. 


P (n ' |n) 


o  (n  -*•  n 1 ) 

T.  o(n  -+  n') 
n' 


(7) 


If  there  were  only  one  independent  observable,  say  the  average  of  the  variable 
x,  then  the  Maximum  Entropy  Principle  would  predict  the  form^’^ 

P(n' |n)  =  P°(n' |n)  exp(-AQ  -  Ax)  ,  (8) 


where  A  and  A  are  Lagrange  multipliers.  P°(n’ln)  may  be  obtained  via  Eq .  (7) 
with  o  replaced  by  a°. 

For  the  present,  we  assume  that  in  an  inelastic  collision  process  there  is 
only  one  independent  observable  which  is  the  average  absolute  fractional  energy 
transferred.  That  is, 


x  =  f  ,  -  f 
n  n 


Af 


with 


f  =  E  /E, 
n  n 


(<>) 


and  therefore, 


c(n  n' )  =  A 


,,  ,,  -A | Af | 

(n  n  )e  1 


2  0°(n  ■+  n'  ) 

n' 


(  U;'* 


where  A  =  o(n)  exp(-AQ).  Of  tlie  two  unknown  constants  in  Eq .  (10),  A  and  ■, 
typically  A  is  a  slowly  varying  function  of  energy  and  the  majority  ol  the 
energy  dependence  of  a  is  contained  in  A.  We  assume  A  to  be  independent  of 


energy.  The  parameter  A  will  be  determined  by  a  procedure  known  as  the  Su_m  _KuJ  e .' 

The  concept  of  the  Sum  Rule  is  as  follows:  The  average  fractional  energy 
transferred  from  an  initial  state  n  is  given  by 


(7) 


5 


i 


r 


<Af  > 


(ID 


2  Af  P(n* |n)  . 


Now  when  n  is  the  lowest  state,  E  ,  >  E  for  all  n'  and  thut  A,f>  >  0.  Here 

n'  n 

the  system  tends  to  gain  energy  per  collision.  When  n  is  the  highest  state 

allowed  for  the  given  E,  E^,  <  for  all  n'  and  thus  <Af><0  and  the  system 

tends  to  lose  energy  per  collision.  Clearly,  for  some  n,  say  n  (perhaps  not 

integral),  <Af>  =  0.  That  is,  when  the  initial  energy  is  E^,  ,  the  average  . 

n 

energy  transferred  is  zero,  the  system  neither  tends  to  gain  nor  lose  energy 
per  collision.  Following  the  lead  of  Procaccia  and  Levine we  assume  that 
E  ~  is  the  microcannonical  equilibrium  value  corresponding  to  an  average  over 
E  h,  namely 


[n  E  1 
L  n  n 


f(n|E) 


(12a) 


where 


f (n I E)  =  (E  -  E  )  g  /Z(E  -  E  )g  . 


(12b) 


Therefore,  using  Eqs.  (4-12),  A(E)  is  obtained  from  the  solution  of 


0  =  ZAf  g  ,(E  -  E  ,)1/2  e~AlAf  I 


When  there  are  two  degrees  of  freedom,  say  and  n,,  then  the  probability 
of  energy  transfer  from  njn2  to  ni'n2'  may  written  as 


P(n! * n2  *  I n j n2 )  =  P(n1'|n1n2)  P(n2' |nj 'n^) 


•?#  ''H ■  v 


where  the  second  term  is  the  probability  of  ending  in  n 2  given  a  particular 
n^ 1 .  The  first  term  is  simply 

P(n1'|n1n2)  =  2  PO^ 'n2*  |  n^)  .  (15) 

n2 ' 


The  energy  transfer  variables  corresponding  to  this  division  are 


Af, 


V 


(16a; 


and 


Af 


2 


(i6b; 


Thus  Af^  is  the  fractional  energy  transferred  to  the  n^  degree  of  freedom  and 
Af2  is  the  fractional  energy  transferred  to  the  n2  degree  of  freedom  given  the 
n^  initial  and  final  states. 

The  state-to-state  cross  section  is  then  given  by 


o(n  n  -+  n  *n  ') 


0(n.n.)  P(n  *n  ' |n  n  ) 


Corresponding  Sura  Rules  are  given  by 


0  = 


and 


2  AfjP^  ' 


1  i 

i/o  1  I Af , 

2  2  Af  g  ,g  ,(E  -  E  -  E  ,yU  e 
n;  n^  1  nl  n2  nl  n2 


(18) 


0  = 


2  P(n1'|n.n2>  £  Af2  P(n2' |  nj 'n^) 

n  '  n0' 


2  2  Af  g  ,g  ,(E  -  E  -  E  ,)  e 

z  nl  n9  ni  nO 

n  *  n  •  1  1  1  2 

1  2 


1/2  -  ‘2|if2l 


.  (19) 


In  Eqs.  (18)  and  (19), 


Af, 


nr 


nl 

IT 


and 


Af 


2 


(20a) 


(20b) 


The  generalization  to  an  arbitrary  number  of  degrees  of  freedom  is  straight¬ 
forward  and  will  not  be  presented  here. 


h 


3.  EXCITATION  CROSS  SECTIONS 


The  only  information  needed  to  apply  the  method  of  Section  2  to  a 
particular  system  is  a  model  for  the  energy  levels.  For  the  present  treat¬ 
ment  we  will  uncouple  all  vibrational  modes.  Moreover,  due  to  the  compute 
time  involved  in  molecules  such  as  CO2  and  H^O  with  many  internal  degrees 
freedom,  we  will  assume  that  hot  bands  do  not  play  a  role  in  the  resulting 
spectra.  That  is,  the  v  and  modes  are  adiabatic  with  respect  to 

excitation.  As  will  be  seen  at  least  for  CO^ ,  this  is  not  a  severe  limitai 
It  should  be  noted,  however,  that  the  present  treatment  may  be  extended  to 
include  all  degrees  of  freedom.  Our  mode]  is  in  effect  placing  additional 
constraints  on  the  system. 

The  energy  level  scheme  adopted  for  CO^  is 


=  w3(v  +  'j)  +  w3x33(v  +  'J2  , 

Ej  -  Be  j(J  +  1)  , 


and  for  1^0 

? 

E  =  oj  (v  +  h)  +  u>3x33(v  +  '5)  ~ 
V3 


E. 

J 


+  ce)  j(j  +  1)  . 


The  spectral  constants  used  are  given  in  Table  I. 
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t  i  ona  1 
dI 


(21a) 

(21b) 

(22a) 

<22h) 


TABLE  1  -  SPECTRAL  COEFFICIENTS^  FOR  C02  AND  H20 


-  — -  , 

co2 

h2o 

W3 

2396.4 

3935.6 

X33 

-  12.63 

-  46.37 

B 

e 

0.39163 

14.575 

C 

e 

— 

9.499 

(a) 

Energies 

in  Units  of  cm 

Employing  Eqs.  (18)  and  (19),  and  X^  may  be  obtained  for  CO 2  and  H^O. 
These  values  are  given  in  Figs.  1  and  2.  Total  cross  sections  are  obtained 
by  summing  over  final  rotational  states,  i.e.. 


o  (0, 

0, 

0, 

0  ->  0, 

0, 

S  0 

(0, 

0, 

0,  0  •> 

o,  o,  V  j’) 

(23) 

r 

A  total  photon  yield  cross  section  may  be  obtained  by 

o  =  2  V'^V)  (2A) 

V 
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AL-80-J7S 
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2  4  6  8  ]0 

RELATIVE  VELOCITY  (KM/SEC) 


Figure  1.  Plot  of  the  C02  Surprisal  Parameters 
(Lagrange  Multipliers)  vs  Relative 
Velocity. 
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This  quantity  is  of  prime  interest,  in  the  plume  signature  since  at  high 
altitudes  the  meanfree  path  is  long  enough  to  ensure  complete  radiative 
relaxat ion. 


Figures  3  and  4  show  plots  of  the  calculated  energy  dependent  cross  sec¬ 
tions.  Solid  and  dashed  lines  show  an  arbitrary  smoothing  in  regions  which 
seem  unphysical.  The  magnitude  of  the  present  cross  sections  were  chosen  to 
agree  with  experimental  values  at  relative  velocities  shown  by  vertical 

arrow.  Unfortunately,  H.,0  excitation  cross  sections  have  only  been  obtained 
at  one  relative  velocity  eliminating  further  comparison  at  this  time.  The 
curves  labeled  "PREVIOUS  a  "  represent  cross  sections  obtained  by  the  author 


using  a  less  exact  treatment 
by  the  present  cross  sections 


(8) 


These  previous  results  have  been  superseded 


Experimental  CO.,  excitation  cross  sections  are  better  known.  Subbarao 

et  al.^^  have  measured  what  are  probably  (0  ■*  1)  ahso  !'..!■•  excitation  cross 

(2) 

sections  and  Rahbee  et  al.  have  measured  relative  (0  >  1)  excitation  cress 

sections.  Figure  5  shows  a  Log  O  vs.  E  comparison  of  the  present  calculation 

with  the  Subbarao  et  al .  data  (points).  Open  circles  indicate  less  accurate 

data.  All  experimental  data  is  shown  for  comparison  on  a  linear  plot  in 

Fig.  6.  It  is  noted  that  a  seems  to  be  a  better  fit  to  the  data.  However, 

1> 

the  experimental  configuration  is  such  that  at  most  one  photon  per  excitation 
collision  may  be  detected. 

In  conclusion,  it  should  be  noted  that  the  theory  presented  here  extends 
the  current  literature  and  thus  more  research  is  necessary  (both  experimental 
and  theoretical)  to  understand  any  restrictions  which  must  he  considered. 
However,  these  results  are  certainly  encouraging. 


(8) 

'Gersh,  M.F.,  Elgin,  .1.,  F.aist,  M.,  and  Bernstein,  I..S.,  "High  Al  tit  tide 
Rocket  Plume  Radiation  Calculations,"  Presented  at  AFCl.  High  Altitude 
Plume  Workshop,  Han scorn  AFB,  MA  (1978). 
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Figure  4.  Calculated  Energy  Dependent  Cross  Sections  for  H.,0 


C02^V3^  T_V  C0LLISI0N  CROSS  SECTION  (CM  ) 


r 


RELATIVE  COLLISION  ENERGY  (eV) 


Figure  5.  Comparison  of  Present  Calculation  with  Data 
of  Subbarnn  et  al P ) 
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